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Abstract 

A class of Fourier Integral Operators which converge to the unitary 
group of the Schrodinger equation in semiclassical limit e — » is con- 
structed. The convergence is in the uniform operator norm and allows 
for an error bound of order 0(e 1 ~ p ) for Ehrenfest timescales, where p can 
be made arbitrary small. For the shorter times of order O(l), the error 
can be improved to arbitrary order in e. In the chemical literature the 
approximation is known as the Herman-Kluk propagator. 



1 Introduction 

We study approximate solutions of the semiclassical time-dependent Schrodinger 
equation 

iej f r(t) = ~Ar(t) + V(x)r(t), ^(0) = VoeL 2 (M d ,C) (1) 

in the semiclassical limit e — » 0. The operator H e := — ^-A + V(x) on the 
right-hand side of |T]) is the so-called Hamiltonian, a self-adjoint operator on 
L 2 (R d ,C). It is well- known that the solution of (JTJ) can be written as 

where the group of unitary operators e~i H 1 is denned by the spectral theorem. 
The semiclassical parameter e may be thought of as the quantum of action h, but 
there are also situations, where e has a different meaning. One example is pro- 
vided by Born-Oppcnhcimcr molecular dynamics, where equation JT]) describes 
the semiclassical motion of the nuclei of a molecule in the case of well-separated 
electronic energy surfaces and e is the square root of the ratio of the electronic 
mass and the average nuclear mass. In this case, the s in front of the time- 
derivative in fl} is due to a rescaling of time t — t/e. This particular choice, 
the so-called "distinguished limit" (see C068J) produces the most interesting 
results in the semiclassical limit e — > 0. 

To formulate our main result, we introduce the following class of Fourier Integral 
Operators (FIOs): 

I^k'-^Mx) := 1 / e^ t ^-y^u(x,y,q,pMy)dqdpdy 7 (2) 

( 27re )3d/2 J R3d 
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where 



^iliP) = (x K (q,p),z> K {q,p)^j is a C 1 -family of canonical transforma- 



tions of the classical phase space T*M. d = R d x 
• S K (q,p) is the associated classical action 



S K (q, P ) 



dt 



X KT (q,p)-E K (q,p)-(hoK T )(q,p) 



dr. 



• the complex-valued phase function is given by 

(x, y, q, p) = S K± (q, p) + S K ' (q, p) ■ ( x - X^ (q,p))-p-(y- q) 



X-X K (q,p) 



(3) 



• and the symbol u is a smooth complex- valued function which is bounded 
with all its derivatives. 

For this class of operators, the authors previously established an L 2 -boundedness 
result, see |RoSw07] . The central result of this paper reads 
Theorem. Let e~~ H 1 be the propagator defined by the time- dependent Schro- 
dinger equation ([1]) on the time-interval [— T, T] with subquadratic potential V £ 
C°°(R d ,R), i.e. sup xeRd \d%V(x)\ < oo for all aeN d with \a\ > 2. 
Then 



sup 

te[-T,T] 



'iH e t 



< C(T)e, 



where = (X K , ) and u are uniquely given as 

• the flow associated to the classical Hamiltonian h(x,(,) = ^|£| 2 + V(x) 



-X K (q,p)=E K (q,p) 
j t ~*\q,p) = -VV (x K \q,p) 



X K (q,p)=q 
Z K "(q,p) =p 



the solution of the Cauchy-problem 

j t u(t,q, P ) = ~u(t,g,p)tr (z- 1 jZ (>' \q,p) 



u{0,q lP ) = 2 d l 2 . 
The <C dxd -valued function 

z(F K \q,pj) = (iid id)F Kt ( 9) p)t 



—i id 
id 



X? (q,p) - iX; (q,p) + iE K q (q,p) + 5" (q,p), 
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depends on elements of the transposed Jacobian 



F K \q,p) f 



( 



xf{q lP ) SJ*(g,p) 
Xf(q,p) 3£*(?,p) 



o//t* wit/i respect to (q,p). 

The equation for it is easily solved. Its solution is the so-called Hcrman-Kluk 
prefactor 



where the branch of the square root is chosen by continuity in time starting 
from t = 0. We presented a simplified version of our main result. Theorem [5] 
will essentially add three central aspects. First, we will state it for more general 
Hamilton operators, namely certain Weyl-quantised pscudodifferential opera- 
tors. Second, for the Ehrenfest-timescale T{e) = Ct log(e _1 ) the result still 
holds with a slightly weaker bound. Third, the error estimate can be improved 
to e N , where N is arbitrary large by adding a correction of the form Y]^_} e n u n 
to u. As u, the u n are solutions of explicitly solvable Cauchy-problems. 
Whereas there is an abundant number of works on Fourier Integral Operators in 
the mathematical literature, only few of them discuss the relation between FIOs 
and the time-dependent Schrodinger-equation. The first works which apply 
FIOs with real- valued phase function to this problem arc [KiKuSI] and |Ki82| . 
In this case one has to deal with the boundary value problem 



To get uniqueness for its solution one has either to restrict to short times t or 
to impose very strong restrictions on the potential. The same problems are met 
in |Fu79j . where Fujiwara applies a related class of operators without integral 
in the oscillatory kernel to the Schrodinger equation to justify the time-slicing 
approach for Feynman's path integrals. 

The avoidance of this problem is the major advantage of complex-valued phase 
functions. In the non-scmiclassical setting, Tataru shows in Ta04 that the 
unitary group of time evolution is an FIO with complex-valued phase function 
(different from ((^J). He also establishes that the simpler choice u(t, q,p) = 2 d / 2 
leads to a parametrix for the non-semiclassical Schrodinger equation. 
A class of operators related to ^ is used in the works [LaSiOOj and [Bu02j for 
the construction of approximate solutions of the semiclassical time-dependent 
Schrodinger equation. In their case, the kernel consists of an integral over the 
momentum space in contrast to the phase-space integral in our expression 



Moreover, these works only allow compactly supported symbols, which enforces 
the truncation of the Hamiltonian in momentum. Finally there is the work of 
Bily and Robert [BiRoOl] , which treats the so-called Thawed Gaussian Approx- 
imation discussed below. 

In contrast to the mathematical literature connecting time-dependent Schro- 
dinger equation and Fourier Integral Operators, there is an abundant number 
of papers in chemical journals on this topic. Nevertheless, the focus is mainly 




Given x, y £ R d , find p such that X** (y,p) 



= x. 
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put on three approximations: the "Thawed Gaussian Approximation" (TGA), 
the "Frozen Gaussian Approximation" (FGA) and the Herman-Kluk expres- 
sion. Confusingly, in the chemical literature both TGA and FGA do not only 
refer to specific algorithms but they are also used to describe whole classes of 
approximations. For example, the Herman-Kluk approximation is sometimes 
considered as an FGA, whereas the TGA refers both to the time-evolution of a 
coherent state and a Fourier Integral Operator. We give a short formal discus- 
sion of the most important methods in the rest of this introduction hinting to 
related rigorous results. 

The starting point is the following identity, which holds for ip £ L 2 (W d , C): 
where 

n e I \ _ \ -\x-q\ 2 /2e ip-(x-q)/e (r\ 

9{ q ,p)K x ) ~ ( n£ y/4 e e W 

denotes the coherent state centered at (q,p) in phase space T*R d . Within the 
chemical community, equation is heuristically explained as an "expansion in 
an overcomplete set of Gaussians" , but the equality can be made rigorous with 
help of the FBI-transform, consider [MaQ2j. Applying the unitary group of (fTJ) 
to expression (j4|), one gets the formal equality 



(e-^WVg) (s) = ^ (e-i^% <p) ) (x)( 9 { q ^) dgdp. (6) 

Hence, one expects an approximation to the solution of (fT]) if the following 
approximate expression for the time-evolution of coherent states is used in ([6]) 

(e-i^^g^) (x) « [det (xf (q,p) + ixf(q,p))] ~* (7) 



x efS" 1 (g,p) e -{x-X K (q,p))-e K (q,p){x-X K (g,p))/2e ^ {q,p)-(x-X« (q,p))/e 



with 



Kt (q,p) = -i (zf(q, P )+iZf(q,pf) (xf{q,p)+iX;\ qi p)) 



-l 



In the chemical literature ([7|) was first derived in |He75j . For rigorous mathemat- 
ical results consider |Ha85| . |Ha98| or |CoRo97] . As the coherent state changes 
its width, expression and the resulting operator were baptised "Thawed 
Gaussian Approximation" . 

However, it turns out numerically (see e.g. the computations in [HaRoGr04 ) 
that more accurate approximations are obtained if one drops the time-dependent 
spreading and uses expressions like ([2]). In the simplest case, the symbol u = 1 
is held constant in t, q and p. This approximation is known as the "Frozen 
Gaussian Approximation" and holds only for times of order O(e), see the remark 
after Theorem[5] To get to the longer times of order O(l), the more sophisticated 
choice of u(t,q,p) as the Herman-Kluk prefactor is needed, sec [HcK184 for the 
original work and [Ka94] and |Ka06| for works, which are methodically related 
to our presentation. Moreover, the latter of them presents the first derivation 
of the higher order corrections. 
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Organisation of the paper and notation 

The paper is organised in the following way. Section 2 will set the stage for 
the discussion of our approximation. Here we will recall central definitions and 
results on Fourier Integral Operators, first and foremost their definition and 
well-definedness on the functions of Schwartz class as well as their bound as 
operators acting on L 2 (R d ,C), see Definition [5] and Theorem [TJ Most of the 
results of this section can be found in |RoSw07j and we refer the reader to 
that paper for a more detailed discussion and motivation of them. In Section 3 
we will prove results on the composition of Weyl-quantised pseudo-differential 
operators and Fourier Integral Operators, see Proposition [2j Moreover, we will 
investigate the time-derivative of a C 1 -family of Fourier Integral Operators in 
Proposition [3] These results will lead to our main result, which we will state in 
Theorem [51 

We close this introduction by a short discussion of the notation. Throughout 
this paper, we will use standard multiindex notation. Vectors will always be 
considered as column vectors. The inner product of two vectors a, b £ K d will 
be denoted as a ■ b — Y2j=i a jbj an d extended to vectors a, b E C d by the same 

formula. The transpose of a matrix A will be A\ whereas A* := A 1 denotes the 
adjoint and finally ej will stand for the jth canonical basis vector of R d or C d . 
For a differentiable mapping F E C 1 (M d ,C d ), we will use both (d x F)(x) and 
F x (x) for the transpose of its Jacobian at x, i.e. ((d x F)(x))jk = (F x (x))jk = 
{d X] F k ){x). This leads to the identity d x (F ■ G) = G X F + F X G for F, G E 
C 1 (R d ,C rf ). The Hessian matrix of a mapping F E C 2 {R d ,C) will be denoted 
by Hess x F(x). 

For the sake of better readability of the formulae, we will be somewhat sloppy 
with respect to the distinction between functions and their values. As a cru- 
cial example, we will write (x — X K (q,p))v for the function (x,y,q,p) •— ► [x — 
X K (q,p))v(x,y,q,p). 

When dealing with canonical transformations, we introduce the following nota- 
tions for a complex linear combinations of the components: 

Z"(q,p) := (0*)* X«(q,p) + i (9*)^ Z*(q,p) 
T{q,p) := (6*)* X K (q,p) - i (6*)^ 

We want to point out that Z K (q,p) is not the complex conjugate of Z K (q,p) for 
non-real matrices O^ . The matrix square root of a positive definite matrix will 
always be chosen as the unique positive definite square root. We want to point 
out that both the determinant of this matrix-square root and the square root 
of a determinant will appear in this paper. 
We define z := <d v q + ip, d z := (6 y ) _1 d q - id p and 

d d 

dw z x(q, P ) = £ (e^ 1 d^x^p) -i^d^x^p) 
k=i j=i 

for functions X E C 1 (E 2d ,C d ), regardless whether they are row or column 
vectors. With these definitions the identity div z X(q,p) — trX z (q,p) still holds. 
Finally, we mention that the expression -^X k{ ' 1 (q,p) ■ S K< ' (q,p) denotes the 
inner product of ^X K< *' S> (q,p) and S k14-s) (q,p). 
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2 Canonical Transformations and Fourier Inte- 
gral Operators 

In this section, we specialise the central definitions and results of |RoSw07] to 
the case of Hamiltonian flows. 



2.1 Symbol classes and canonical transformations 

The definition of our FIOs involves two fundamental objects. One of them is a 
smooth complex-valued function, the so-called symbol. The following definition 
deviates from [RoSw07j by the additional e-dependence. 

Definition 1 (Symbol class). Let m = (mj)i<j<j G K.' 7 and d — {dj)i<j<j G 
N J . We say that u :]0, 1] x Rj d — > is a symbol of class S[m; d], if there is 
Eq < 1, such that u e G C°°(IR.' d ', C N ) for all e < Eq and the following quantities 
are finite for any k > 



M™[u] := sup max sup 

e<e l«l= fe 2GRl d l 



J 



(8) 



where (z) := y/l + \z\ 2 . We extend this definition by setting 

5[+oo;d]= |J ... |J S[( mi ,...,mj);d}. (9) 

mi£R mj£l 

The second central object in the definition of a Fourier Integral Operator is a 
canonical transformation of the classical phase space. 

Definition 2 (Canonical transformation). Let n(q,p) = (X K (q,p), E K (q,p)) be 
a diffeomorphism of T*R d — M. d x M. d . We represent its differential by the 
following Jacobian matrix 



F K (q p)=( X « iq > p) ! X P^1) (10) 

k is said to be a canonical transformation if F K (q,p) is symplectic for any 
(q,p) in T*R d , i.e. 

F K {q,p) G Sp(2d) := {S G G\{2d)\SKjS = J) with J^f^ ^ . 

To get good properties for our operators, we need to restrict the class of canon- 
ical transformations under consideration. 
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Definition 3 (Canonical transformation of class B). A canonical transforma- 
tion k of T*R d is said to be of class B if F K € S[Q;2d]. A time- dependent 
family of canonical transformations K will be called of class B in [-T, T] if it 
is pointwise continuously differentiable with respect to time and we have for all 
k>0 



sup Ml 

te[-T,T] 



F H 



< oo ano 



sup M" 

te[-T,T] 



dt 



F H 



< oo. 



In particular F K and -^F K are of class S[0;2d] pointwise for t 6 [— T, T]. 
We also have to restrict the Hamiltonians we use. 

,C°°(R 2d ,C)) is called 



Definition 4. A time- dependent Hamiltonian h G 
subquadratic, if 



sup sup \\d? ft(*, at, 

-T<t<T (x,J)6R <i xR d 



(11) 



is finite for all \a\ > 2 and T > 0. is called sublinear, if the quantity is finite 
for all \a\>l. 

The next result will investigate the relation between classical Hamiltonians and 
the flows they generate. 

Proposition 1. If h € C(R, C°°(R 2<i , C)) is a time- dependent subquadratic Ha- 
miltonian, the Hamiltonian flow k^'^ generated by h, 



j t K^ = JV (x , S) M*,k (M) ), k (s ' s) = id 



(12) 



is a family of canonical transformations of class B in [-T, T] . Moreover, every 
Hamiltonian flow of class B is generated by a subquadratic Hamiltonian. 
Under the additional assumption £)ft>||£°o(KxK 2,i ) < oo for all 2 < \a 
no + 2, we have 



< 



sup Mj: 

\t-s\<T(e) 



<C k (2C T ) k \\oge\ k e- 2K ^, 



for all k < no on the Ehrenfest timescale T{e) — Ct loge , where 
Kk(T) = ma,x sup sup \\dt 6) }iess {x4) h(t, x, 

|a|=ft _ T <t<T ( a; ,£)eR'*xK d 



Proof. The basic identity follows by differentiating (fT2")l with respect to (g, p) : 
^ W> (q,p) = me SS(x ^h(t, K^(q,p))F^ 3) (q,p). (13) 
The fundamental theorem of calculus gives 



F Kit - 3 \q, P ) 



< 



Ress (x!i) h(T,K^ s \q, P )) F^'* ] (q,p) 



dT 



lidl 



<K*(T) 



F h 



7 



Hence, the bound for the Jacobian follows from Gronwall's Lemma: 

< e K%{T)\t-s\_ 



For the derivatives we have 

t 

< 



d? p) F^ s \q,p) 



K»{T) d? p) F^"{q,p) 



dr 



E 

f3<a 



%,P) 



Ke SS{x , h(T,K^ s \q,p)) \ d^ p) F K ^'(q,p) 



dr 



so Gronwall's Lemma provides inductively 

^ )P) F« Ct ' s, (g,p)|<C fc (2T) fe e^mi*-l 

where Ck depends on Kf(T) for I < k. The result for the Ehrenfest timescale 
follows by substituting T(e) — Ct log(e _1 ) into this expression. 
Now consider a Hamiltonian flow of class B. The identity (fT5|) gives 

.(*,') t 



J ( -F KK ' '(q,p) ) J | F ' '(q,p)) J = EeBB M h{t,KS t <'>(q,p)) 



Hence, h is subquadratic, as ^jF k( is of class S[0; 2d] by definition. 
Remark 1. 



□ 



1. It is easy to show that there is an equivalence between linear Hamiltonian 
flows and quadratic Hamiltonians. Every quadratic Hamiltonian generates 
a linear flow and every Hamiltonian linear flow K^ t,s ^(q,p) = M(t, s)(^), 
M(-, s) £ C 1 (IR, Sp(2(i)) is generated by a quadratic Hamiltonian, namely 

h M (t,q,p) = l -(q p) (j (j^M(M)) JM*(t,s)j) Q. 



2. By estimating the logarithm, we can have a bound of the form 



sup Mfc 

l*-*l<T(e) 



(14) 



for the Ehrenfest timescale, where p(Ct) < Po for any Po > 0, if Ct is 
chosen small enough. 

3. From now on, all considered canonical transformations are assumed to be 
of class B. 

An important quantity associated with a canonical transformation is the so- 
called action. 

Definition 5 (Action). Let n(q,p) = (X K (q,p), S K (g,p)) be a canonical trans- 
formation ofT*M. d . A real-valued function S K is called an action associated 
to k if it fulfills 



S*(q 7 p) = -p + X«(q,p)Z K (q 7 p), S«(q,p) = X*(q 7 p)Z K (q,p). (15) 
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Remark 2. 

1. An action associated to a canonical transformation is only denned up to 
an additive constant. If we consider a time-dependent family of canonical 
transformations k* , we will choose this time-dependent constant such that 
S K (liP) is C 1 with respect to time. 

2. If k(*' s ' is induced by a Hamiltonian h(t,x,£), the action of classical me- 
chanics 

t 

Sf s) (q,p)= J (±X" T (q,p) - E« T (q,p) - h^K^fap))} dr 

s 

is an action in the sense of this definition. In this case, we use the conven- 
tion S K ' (q,p) — 0, where 5 K< ' (q,p) is now considered as a function of 
t. We cannot assume S ld (q,p) = 0, as the case h(t,x, £) = h(t) shows. 



2.2 Definition of FIOs and continuity results 

In this section, we will define the operators we will use to approximate of the 
propagator of the Schrodinger equation. 

Definition 6 (Fourier Integral Operator). For u £ 5[(-foo, m p ); (3d, d)}, a 
Schwartz-class function ip £ 5(R d , C) and n > vnP + d, we define 



[l s { K ;u;Q x ,@y)ip]{x):= (16) 
-J— f ei^^y^ e ^ eV HLl) n [u(x,y,q,pMy)]dqdpdy, 

where 

• Q x and Q v are complex symmetric matrices (i.e. = with positive 
definite real part, 

• the complex-valued phase- function is given by 

$ K (x, y, q,p; Q x , 0*) = S K (q,p) - p ■ (y - q) + E K (q,p) ■ (x - X K (q,p)) 

+ ^(y-q)-Q y (y-q) 

+ ^(x-X K (q,p))-Q*(x-X K (q,p)) 

• and the differential operator L y is defined by 
1 



Li, 



1 - ieV y ^(x, y, q,p; Q x , Qv) ■ V y 



y l+\V y $K(x,y,q,p;Q*,Qy)\ 2 
The operator L y , the Banach-space adjoint of L y , i.e. 

v(y)[Llu](y) dy = [L y v}(y)u{y) dy, 



makes the integrals absolutely convergent. Moreover, it is constructed such that 
expression (|16j) coincides for different n with n > m p + d (especially with n = 0, 
if m p < -d). 

The following theorem combines the central results of |RoSw07] . 
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Theorem 1. 

1. If u G 5[+oo;4d], T e (k; u; Q x , Q v ) sends <S(R d ,C) into itself and is con- 
tinuous. 

2. If u G S*[0;4(i] 7 T e (k;u]@ x ,O v ) can be extended in a unique way to a 
linear bounded operator L 2 (M. d , C) — > L 2 (R d , C) and f/iere exists a constant 
C(Mg,Q x ,Q v ) such that 

\\F{K-u-Q\Qy)\\ L2 ^ L2 <C{M£-Q\Qy) ll^x,,)«[U- (17) 

|a|<4d+l 

In i/ie special case where u G SfO; 2d] is independent of (x,y), we have 
\\T{n;u-Q x ,Q y )\\ L2 ^ L2 < 2- d/2 dot {?Ste x m y y Li \\u\\ L aa. (18) 

Remark 3. 

1. The dependence of C(Mq ; G 31 , & y ) on Mq , 0^ and 6 y can be made more 
explicit. Consider [RoSw07j for the precise expression. 

2. There is an analogous result for Weyl-quantised pseudodifferential opera- 
tors 

(of(hW)(x) := l^'^h (^T^) ^dydi, 

see for example (Ma02 : 

(a) If h G 5[+oo; 2d], op £ (h) sends <S(M d , C) into itself and is continuous. 

(b) If h G S[0; 2d], op e (h) can be extended in a unique way to a linear 
bounded operator L 2 (R d ,C) -> L 2 (R d ,C) with e-independent norm 

\a\<2d+l 

The second part is the famous Calderon-Vaillancourt Theorem. 

3. We have T E ^id; det(Q x + O a )5; O 2 , ®y\ = id, compare the appendix for 
the correct choice of the square root. 



3 Composition with PDOs and time-derivatives 

The standard approach in the field of asymptotic analysis consists in a two step 
procedure. First, one constructs an asymptotic solution Uff(t,s)ipg of order 
0(e N+1 ), i.e. a function which fulfills 

(iej t -H E (t)^ Uh(t,sW s =e N+1 R%(t,sW s . (19) 

If one can establish an e-independent bound on the remainder R E N (t,s), the 
asymptotic solution can be turned into an approximate solution of the unitary 
group with help of a special version of Gronwall's Lemma, (see for example 
Lemma 2.8 in |Ha98] for the strategy of the proof): 
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Lemma 1. Let U e (t,s) be the propagator of the time- dependent Schrodinger- 
equation 

(iej t - H*(t)\ r(t) = 0, r(s) = r s G D c L 2 (R d , C) 

for some family of self-adjoint operators H e (t) with common domain D. More- 
over, for some T > and —T < t, s < T let U^{t,s) be a family of bounded 
operators, which is strongly differ entiable with respect to t, leaves the domain of 
H e {t) invariant and which fulfills 

izjUHt, ~ H e (t)U £ N (t, s)r(s) = e N+1 R%(t, s)ip e (s) 

with Ufj(s, s) = id. If \\R E N (t, s)||i2_ >i 2 < oo for all —T < t, s < T, we have 

\\U e N (t,s)-U e (t,s)\\ L2 ^ L2 <s 

In this section, we state the intermediate results needed for the construction of 
the asymptotic solution. 

In Proposition we show using Weyl-quantisation that the composition of 
differential operators with Fourier Integral Operators is again an FIO. Moreover, 
we give an asymptotic expansion of the symbol of the new FIO, whose terms 
but for the last are x- independent. This is important, as x-dependence of the 
symbol may be converted to e-dependence, which can be seen from Lemma [3] 
Proposition [3] deals with the time-derivative of a family of FIOs. Finally, we 
will establish an uniqueness result for symbols and canonical transformations in 
Proposition [H 



t 

J \\R e N (T,s)\\ L ^ L 2 dr . 

s 



3.1 Statement of intermediate results 

i_ 

To state our results, we need the matrix Z(q,p) — Z^(q,p) (Q x ) 2 , which already 
appeared as Z (F K (q,p)) in the statement of our main result in the introduction. 
We justify this abuse of notation by better readability of the formulae presented 
here. The invertibility of Z(q,p), which is implicitly claimed in the following 
statements, is shown in Lemma [2] 
The composition result reads: 

Proposition 2. Let h £ S[mh, 2d] be polynomial in £ and u £ S[m u ; 2d]. Then 
we have 

op s (h)I s (K;u;e x ,@y)=l £ ^K;J2s n v n ;e x ,Q^j + e N+1 T (k; v £ n+1 ; 9 x , 0») 
as operators from S(R d ,C) to S(R d ,C). 

v n £ S[m u + mh', 2d], n < N and v e N+l £ S[(rn,h, m u + m^); (d, 2d)] are given by 

v n (q,P) = L n [h; k; Q x , Q v ]u(q,p) 
v £ N +i(x, q,P) = L s N+1 [h; k; Q x , e v ]u{q,p), 
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where L n [h; k; Q x , Q v ] and L e N+1 [h; k; Q x , Q v ] are linear differential operators 
of order n and N + 1 in {q,p). The coefficients of the L n [h; k; & x , <d y ] are 
rational functions, with a numerator depending on derivatives of h from order n 
to 2n and derivatives of F K (q,p) of order < n, and a denominator of the form 
det(Z(q,p)) m for some m > 0. L £ N [h; k; @ x , O v ] is of the same form depending 
on derivatives of h from order N to 2N + 1 and derivatives of F K (q,p) of order 
<2N+1. 

The explicit expressions for v ,vi and v 2 are 
v (q,p) =u(q,p)(hoK)(q,p) (20) 

Vl (q,p) = -div z [{(h x + iQ x h s ) o K (q,p)) 1 Z-\q,p)u{q,pj) (21) 
+ u(q,p)^tl (Z- 1 (q,p)d z ((h x + iO x hs) o K(q,p))) 

v 2 (q,p)=L 2 [h> 3 ;K;e x ,e y }u(q,p) (22) 
1 d 

+ 2 E div * ( u ^P) d ** + ^ X difh) ° K(q,p)Z-\q,p)e k )] f Z-\q,p)) 
k=i 

The linear partial differential operator L 2 {h> 3 ; k; Q x , Q y ] depends on derivatives 
of h of order 3 and 4. 

Remark 4. As the coefficients of the differential operators are rational functions 
of the form described in the statement, a bound for the elements of F K (q,p) 
and their derivatives of the form Ce~ p gives a bound for the coefficients of the 
form C's~ Mp for some MeN, where C depends on derivatives of h. 

The second result of this section will investigate the time-derivative of a family 
of FIOs. In the case of a time-dependent family of canonical transformations, 
we have the following result: 

Proposition 3. Let u G C(R, S[(m q , m p ); (rf, d)]) be a family of time- dependent 
symbols with u(-,q,p) G C 1 (R, C) and (-^u)(t, •,•) G S[(m q , m p ); (d, d)], k* a 
family of canonical transformations of class B, S K an action associated to k' ; 
Q x G C 1 (R, Gl(c?)) a family of complex symmetric matrices with positive definite 
real part and Q v complex symmetric with positive definite real part. We have 

ie^-F («*; «; B x (t), 0") = I £ ( «*; £ s n v n ; Q x (t), qA 



n=0 
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with 



v (t,q,p) = u(t,q,p) y- 
vi(t,q,p) = i^-u(t,q,p) 




(23) 



(24) 




) 



v 2 (t,q,p) = -Y,div z [d Zk (-e x {t)Z- 1 (t,q,p)e k u(q,p)j Z-\t,q,p) 



where vq, fi, «2 G C° (R, <S[(m q , m p ); (<i, d)]). 

Remark 5. In both propositions, the case of a linear canonical transformation, 
a quadratic symbol h and a constant symbol u results in v n = for n > 2. 
This will result in the exactness of the Herman-Kluk expression for quadratic 
Hamiltonians. 

Finally, we have the following uniqueness result. 

Proposition 4. Let K\ and k 2 be two canonical transformations of class B and 



then u = v and Ki(q,p) = K2{q,p) for all (q,p) E suppu. 

3.2 Some auxiliary results 

The proofs of these results will strongly rely on results about conversion of 
x-dependence to ^-dependence. We start with the following auxiliary result. 

Lemma 2. We have 



i$ K z (x, y, q,p; B x , 0") = Z*{q,p) (0*)* (x - X K (q,p)) =: Z(q,p)(x - X K (q,p)). 



Z(q,p) = (i(<d y ) 1 id)(F K (q,p)y(— iO x id)t is invertible and its inverse 
Z^ 1 (q,p) is in the class S[0;2d\. 



(25) 



u,v e S[0;2d\. If 



lim \\r( Kl ;u;e x 7 e y )-i e {K 2 ;v : e x ,e y )\\ 

£^0 



= 
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Remark 6. 

1. We recall that 

z K z { q , P ) = (Vr 1 d q - id p ) ((e*)' + * (e*r* s K (<z, P )) . 

2. Obviously, Z(q,p) depends on q and p only via the elements of F K (q,p). 
For better readability, we do not explicitly denote this dependence. More- 
over, we drop the dependence on <d x and Q v in the notation. 

3. For a linear canonical transformation n(q,p) = M , with M 6 Sp(2d), 

we have F K (q,p) = M, so Z(g,p) = (iid (0")- 1 )Mt(-ie ai id)+ is 
constant with respect to (q,p). 

Proof. The derivatives of $ K (x, j/, g,p; 8 X , O y ) with respect to q and p are 

y, q, P] e x , e») = [Hg - l x g K e*](g,p)(z - - *e^(y - <?) 

y, q,p; G x , 0") = [S£ - tX£e*](«,p)(a: - X K (g,p)) - (y - g), 

which gives the identity for Z(q,p). Obviously, Z{q,p) inherits its symbol class 
from F K (q,p). Moreover, we have 

Z(q,p) {®& X )~ X Z(q,p)* 
= 25R(9^)- 1 + (A(Q x )F K (q,p) ( i(9 ^ ')) U(e x )F K (q )P ) (* * 
with 

Hence, by the superadditivity of the determinant for positive definite hermitian 
matrices, det Z(q,p) is uniformly bounded away from for all q and p, so by its 
expression via the formula of minors, Z~ 1 (q,p) G S[0; 2d], as Z(q,p) is. □ 

We introduce the following notation: 

Definition 7. Two symbols u,v G 5[+oo;4d] are said to be equivalent with 
respect to k if 

u; Q x , Q v ) = T(k; v; Q x , Q v ) 
as operators from S(M. d , C) to S($L d , C). In this case we write u ~ v. 

The central technical identity, on which most of the following results rely, is con- 
tained in the following lemma, which is a special case of Lemma 5 in |RoSw07] . 
Here we present a proof based on definition whereas in |RoSw07| an alterna- 
tive definition based on a smoothing of the oscillatory integrals is used. 

Lemma 3. Let u S S[(m x , m q , m p ); (d, d, d)] . Then 

(xj - Xf(q,p))u ~ ev, 
where v 6 S[(m x , m q , m p ); (d, d, d)] is given by 

v{x,q,p) = -div 2 ^ J Z^ 1 (q,p)u(x,q 1 p)^ 



14 



Proof. We introduce a countable, locally finite partition of unity xi{QtP) with 



C < oo. Let m > m p + d. Using dominated 



convergence, we get 
J £ ( K ; ( Xj -X«(q,p))u;e x ,ey) <p(x) 

1 00 f 

= ~R — / Xl(q,P)(xj - X^{q,p))e^ K '(Ll) m (wp)(x,y,q,p) dqdpdy 

1 °° f 

= 7^ — \3dJ2Yl / Xl(q,P)(xj -X?(q,p))ei* (u<p)(x,y,q,p) dqdpdy 

00 „ d 

= j2^^^J^ d XM^P)^Y^d Zk (Zj~k u ) ( x >Q,P)p(y) dqdpdy (26) 



fc=i 



00 d 

X) / ^2( d z k Xi)(l,P)e i ' S '' i (Z~ k 1 u(p){x,y 1 q 1 p) dqdpdy (27) 



After introducing m powers of ij, , the integrand in (|27p is dominated by 

d 
k=l 

Thus (|27p vanishes, whereas (j2"B"|) is the expected quantity. 



□ 



By iterative applications, the previous lemma easily extends to any polynomial 
x dependence. We will only state the result for quadratic polynomials and refer 
the reader to the proof of Proposition [21 where a similar problem is treated, for 
the general case. 

Proposition 5. Let P(q,p, x) = a(q,p) + a(q,p) ■ x + x ■ A(q,p)x be a quadratic 
polynomial with coefficients in S[mp;2d] and u 6 S[m u ;2d\. Then 

2 

P(q,p,x-X K (q,p))u(q,p) ~ ^V l u„(<?,p), 

n=0 

where the v n G S[mp + m u ; 2d] are given by 
v Q (q,p) = a(q,p)u(q,p) 

vi(q,p) = -div z (a(q,p) Jt Z^ 1 (q,p)u(q,p)) 

+ tT{Z- 1 (q, P )X^(q,p)A(q,p)) u{q,p) 



V-2 



(liP) = X diVz ( d zkl A (<l,P) z 1 (l,P)eku{q,p)]' ! Z 1 {q,p)) 



k=l 



Proof. The statement follows from two successive applications of Lemma[3] The 
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quadratic part contributes to v\ and Vi\ 

(x - X K (q,p)) ■ A(q,p) (x - X K (q,p)) u(q,p) 
~ -ediv z (u(q,p)(x - X K (q,p))^A(q,p)Z- 1 (q,p)) 
= -etr (d, [u(q,p)(x-X^q,p)^A(q,p)Z- 1 (q,p)]) 
= etr (A(q,p)Z- 1 {q,p)X*(q,p)) u(q,p) 



-e^x - X K [q,p))^ (^(AfapjZ- 1 (9,p)e fcU ( 9 ,p))) 
fe=i 

~ etrtZ-^g.pJXJfapJAfop)) u(g,p) 

d 

+ e 2 ]Tdiv z [d Zk (A(q 7 p)Z- 1 (q 7 p)e k u(q,p)yZ- 1 (q,p)] , 
fc=l 

where we used the invariance of the trace under cyclic permutations. □ 
3.3 Proofs of the Propositions [2HH 

We have now established all auxiliary results necessary for proving our central 
propositions 

Proof, (of Proposition^ Let ip 6 iS(M d ,C). The composition of op £ (h) with the 
FIO applied to ip is 

[op e {h)T{K]u;® x ,® v )ip]{x) 



(2^)5^/2 J RSd 

where ^> K (x, w,y,q,p; <d x , <d v ) := £ ■ (x — to) + <5> K (w, y, q,p; <d x ,Q y ). 
We introduce the following combinations of positions and momenta, which cor- 
respond to creation and annihilation "variables" 

a(x,0 := {&>)$x + i(&°)-$t, 
a(x,£) := (e*)iz-i(e*)-i£ 
and their "dual operators" 

d a :=l((G*)-id x -i 

The Taylor-expansion of the symbol h to order 2N around n(q,p) reads 
M*>0 = E ^ {{ d a4h) °") (?,P)(a-Z«(?,p)) a (a-^(?,p) 

|a+/3|=2JV+l 



fc, T fa - Z K (q,p),a - Z K (q,pf) + h R (a,a,q,p) 
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where 

Ra.ji (a,a,q,p) 



= / v^®- 1 (dSdgh) {x + v(X«(q,p)-x),Z + v{E«(q,p)-S))d*. 

o 

We split the integral into a part which contains the Taylor polynomial /it and a 
remainder containing hji. In the first step, we discuss only the part containing 
/it- We have, dropping the argument (^4^, £) of a and a, 

i (0*)* % - i (0*)"* d w ) ^ = a - Z K (q,p). 



Moreover 

\ (0*)* d € - i (0*)-* ^ (a - Z"(g,p)) = 

i (0*)* 8 £ - i (0*)"* d w ^j (a - Z K ( g ,p)) = -i id. 

By integration by parts we get 

(a-Z K (q,p)) [a-Z(q,p)j v{q,p) ~ — - — [a - Z {q,p)J v(q,p), 

(28) 



where we extended the meaning of "~" in an obvious way. We have 

Q (®T 2 d ( + i (0 X )"^ d w + 2i (0*)* iT^p)^ (a - Z*(«,p)) 

and hence 

£ div z (eji(e a )^- 1 ( g ,p)(a-^ t ( g ,p)) 7 " e \(g,p)) 
■jpj 2^ [ 2-^(l- e k) m [a-Z (q,p)j (C {ek>em) v){q,p) (29) 



fc|7fc#0 \m=l 
+ 



(a- Z (q,p)J (C ek v)(q,p) J , 
where the differential operators £( efe , em ) and C ek are given by 

(£ {ektem) v)(q,p) := 2el (0*)* Z-\q,p)dX W&p) 
(C ek v)(q,p) := -2div z (e* (0 X )' iT^pM^p)) 
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and #7 denotes the number of non-zero components of 7. 

The symmetrization by the summation over k allows for the iteration of the 
procedure. We define the three sets 

T x := {7 e N d I |7| = 1} , T 2 :=Y x xT x , T^i^UrY 

In expression (|29p the sum is taken over all possible reductions of the multi- 
index 7 by elements of the "brick-sets" Y\ and T 2 - After another integration 
by parts in all terms with (li — Z K (q,p)^j -dependence, the sum is taken over all 
possible reductions of 7 by elements in T x T, which may be considered as a 
two-step path in T, plus the terms which already led to 7 = in the first step. 
So after the removal of all — Z (q,p)j -dependence, the sum is taken over all 

possible paths in the "brick-set" T which reduce 7 to zero. To formalise this 
idea, we define the map 

f ■ 1 : r -> N d 



[7] 

With 

A(7,7l: ■ ■ -)7n) = 



7 7 e Ti 

71 + 72 7 = (71,72) e r 2 . 



#(7- EN)) 7ner! 

Kn J 

#(7~ EN)) (7 - E [7/] - ej) 7n = (e j ,e k )GT 2 

l<n / \ l<7i / u 



we have 



a-Z K (q,p)Y v(q,p) (30) 
~ ^2 e fe A(7,7i, . . . ,7 fc ) . . . A(7,7i,7 2 )A(7,7i) (C lk . . .£ yi v) (q,p). 

[71]+.. .+[7*]=7 

Combining (|2"5j) and (J5U]) , we get 

(h T u)(q,p) 

N 

~ J] e"£„[/i; k; 6 s , e>(?,p) + £ ff+1 I^ +1 [/>; «; 6*, ® y ]u(q,p) 

n=0 

= E E ^3^1 Il A (7>7i,--,70A ( («^*° 

\/3\<2N 7i..., 7fc Gr V=l / 

a</3 [7i]+...+[7fc]=/3-« 

where e Af+1 L^ r+1 [/i; k; O 1 , <d v ]u(q,p) contains all the terms of order e N+1 and 
higher. As k ranges between \\(3 — a|/2] and \/3 — a\, we have 

L n [h;n;Q x ,& y ] 

/n-\a\ 

= E E ^i^.^i II A( 7 ,7i,...,7l)£y, 



alffl - a)! 

n<|a+/3|<2n 7i,---,7„-|"l er ' y i=l 

a</3 [7i]+-+[7»-|«|]=/9-a 
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with the following convention for n — \a\ = 

n— | ct | 

^2 II A (7>7i, ■ ■ ■ ,7i)^7i = id - 

71i'",7n-|a| '=1 

For the first few terms in the expansion, we have more transparent expressions. 
The zeroth order term 

(Hok) (q,p)u(q,p) 

is provided by a = (5 = 0. For the first order term, there are three contributions. 

1. The terms with |/?| = 1, a = /3, which result in 

etr (((d a dah) o k) (q,p))u(q,p). 

2. The terms \j3\ = 1, a = 0, which give 

-ediv, (((^) f o «) ( ? , p ) (0*)* Z- l {q,p)u{q,p)) . 

3. The first order contribution of terms |/3| = 2, a — 0, which is 

etr (z- l {q,p)d z Z K {q,p) ((Hessfch) o n) (q,p) (0*)*) <q,p). 
By an application of the chain rule, they may be combined to 

= -«iiv z (((£U) f ° «) (?,p) (9")* iT^pM^p)) 
+ etr (z^fopjft, ((0*)' ((^A)o«) (g,p))) u(g,p). 
The second order term arises in a similar way. 

The form of the coefficients of the differential operators L n [h; n; <d x , 9 y ] follows, 
if Z^ 1 (q,p) is expressed by the formula of minors. With respect to the symbol 
class of v n , it is sufficient to note that k 6 S[l; 2d]. 

We turn to the discussion of the remainder. The (a — Z K (q,p)) and (a—Z (q,p))- 
factors may be converted to e-dependence analogously to hx, resulting in terms 
of order e N+1 to e 2N+2 . As h{x, £) is polynomial in £, the resulting expression 
equals the application of a differential operator of order to an FIO. We have 
d x & K (x, y, q,p; Q x , Q v ) — S K (<7,p) + iQ x (x — X K (q,p)), hence the symbol class 
by iterative applications of Lemma [31 □ 

Remark 7. In the case k = id, u(q,p) = 1, X = Q v = id, the proof provides the 
asymptotic expansion of the Anti-Wick symbol of a Weyl quantised pseudodif- 
ferential operator. 

We have Z(q,p) = 2 id, Z K = q + ip and Z = q — ip. Hence d z Z = 0, 
£(e-,e h ) — an( i {^-e j u){q,p) = — (d Zj u)(q,p). Moreover, all £ ej commute. 
By straightforward calculation, the Anti-Wick symbol is 

JL (_i)\is\ E m 

J2 { -\^(^)(q,p), 
|/91=0 P - 
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where 

d 

v'_TT, l2 



fe=i 

Thus formally 

^AW = e £A /lWeyl, 

so we recover that the Anti-Wick quantisation is the solution of the Cauchy- 
problem for the inverse heat-equation at time t — e with the Weyl-symbol as 
initial datum (compare Ma()2 . where the transition is expressed by the heat- 
kernel). 

Next, we give the easy proof of Proposition [3] 

Proof, (of Proposition [21) By direct computation is-^X e (k';«) is an FIO with 
symbol 

iej t <t, q,p) - u(t, g,p)^*"' {x, V, q,p) @ x (t), Q v ) 



= i£-jru(t,q,p) - u(t,q,p) 



|^( ff ,p)-|^( ff .p).^fop) 



+ i(x- X*™ (q,p)) ■ |e-(t) (x - X"™ (q,p)) 

The expressions (|23p and (|24p follow from Proposition \5[ □ 

We close this section with the proof of Proposition [U 

Proof, (of Proposition H{ The proof relies on the inner product 

/o £ ',r „Z £ (K]w;<d x ,@ y )g';- e \) (31) 

\ 3 k'(<3o,Po)' v ;i, (9o,Po)/ v ' 

for symbols u> £ S*[0; 2d] and canonical transformations n, k! of class £>, where 



e ,0 ^ = det^G)^ -{x-q)-e(x-q)/2e p ip-(x-q)/s 



Straightforward calculation gives 

_ 2 d det(3?e x )*det(3?e^ f i „ oo1 , , , , 

i n£ ) det(29 a: ) 5 det(2ey) 5 •/ 

where 

^ K '( g ,p,g ,Po) = ^(9,p) 

+ \{V - <7o)(p + Po) - ^{X K {q,p) - X K '{q ,p Q )){E K (q,p) + S K '(g ,p )) 
+ i(q -q)- Q y (q - ?) /4 + i (po - p) • (6T 1 (p - p) /4 
+ * (x K '(q ,p ) - X*{q,p)) ■ e x (x K '(q ,po) - X K (q,p)) /4 

+ i (s K '(g ,p ) - ~ K (g,p)) • (G*)- 1 (s K '( 9 o,Po) - S K (g,p)) /4. 
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We chose a 6 C{j° (M. 2d ,M) with a = 1 in an neighborhood of (qo,po) and split 
the integral into 

(4C,,o)» :r£ ^^ 0X ' 0a )^Co)) ( 32 ) 

+ <.<C,po)^ £ («; (i - © x 5 © a ) ■ 03) 

It is easily seen that 

^{q,P,qo,Po) = and (V( g , p )3?*)(<?,p, ?o,Po) = 



if and only if (<7,p) = {qo,Po) and K,(qo,po) = k'(<7o,Po)- Thus the phase in 
is non-stationary on the support of w(l — a - ), so after integrations by parts with 
the operator 



-is 



\\V( q ,p)^(qo,Po,q,p)\\ 



! V (?,p)*(9o,Po,g,p) • V(,, p ) 



we have lim e ^o (EU) = 0. By the same argument, the case k(<7q,Po) 7^ k'(<70j£>o) 
gives lim e ^o Q32|l = 0. In the case k^OjPo) = ^'(<Zo>.Po)j we have 

Hess ( , 0iPo) ^< K = 5 (e y)-i) + ^ K (9o,Po) f (e-)" 1 ) FK (lo,Po), 

at the stationary points, so by the Stationary Phase Theorem (Theorem 7.7.5. 
in [Ho88] ) 

l^o (^Jpo)' J£ ("! ^ *> ^ 9 E Z, P0 )) = C ^ e "i >(<Zo,Po), 
with the non-vanishing constant 



2d det(Ke a: ) 3 (K9^ 3 



det(9 a; ) 2 det(9y) 2 

_ 1 

det ((7 (e'- 1 ) +FK(9o ' Po)t (T (eVO^*'* )" 

Subsuming this discussion, we have 
0=lim/o £ 'T „[X e (Ki]u;Q x ,Q y )-l £ (K2;v;Q x ,& y )]g e 'f ,) 

_ Jc[ki; 6 s ; 6 a ]u((7o,Po) ki(9o,Po) ^ n 2 (q ,p ) 

\C[ki] Q x ] @ v }u(q ,p Q ) - C[n 2 ; Q x \ ® y ]v(q ,p Q ) Ki{q ,po) = K 2 (q ,p ) 

In the case Ki(qo,po) K 2(<Zo,Po) w e immediately get u(qo,po) = and by 
symmetry v(q ,po) = = u(q ,Po). In the case Ki(q ,p ) = K 2 (g ,Po), we 
either have u(qo,po) = or u(qo,po) ^ 0. In the first case, we immediately 
get v(qo,pa) = = u(qo,po)- In the second case u does not vanish in a 
neighbourhood of (qo,po). Hence n\ — k 2 in the same neighborhood, thus 
C[«i; ® x l ® v ] = C[k 2 - Q x - 9»] and so u(q , Po ) = v(q ,p ). □ 
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4 Uniform approximations of the unitary group 

In this section we will combine the results of Section [3. II to our main result. 

Theorem 2. Let U E (t,s) be the propagator associated to the time- dependent 
Schro dinger- equation 

isj t r(t) = H e (t)r(t), = ft e L 2 (R d ,C) 

on the time-interval —T < s,t < T where H £ {t) = op e (ho + ehi) with sub- 
quadratic ho(t,x,£,) and sublinear hi(t,x,^), both polynomial in £. Moreover 
let 

• Q y £ Gl(d) be complex symmetric with positive definite real part and 

• 6 1 £ C 1 (R, Gl(d)) be complex symmetric fulfilling < 7id < $l® x (t) < 
j' id for all t £ [— T, T] in the sense of quadratic forms. 



Then 



sup 

-T<s,t<T 



U £ (t, s) - T ( «(*•»); e n u n ; ©•(*),©») 

V n=0 J 



< C(T)e 

L 2 ->L 2 



N+1 



where k^' s ^ and the u n are uniquely given as 

• the Hamiltonian flow associated to ho and 

• the solutions of 

d . . 1 . . 
— u n {t,s,q,p) = -u n (t,s,q,p)x 



tr (z-\t, s, Q,P)j t Z(t, *, q,p)j ~ ihi (t, X^ s) (q,p), (q,p)) 
Ediv 2 (d Zk (^jQ x {t)Z-\t,s,q,p)e h u n ^ Z-\t,s,q,p) 



k=l 

n-2 

i 

n— 3 

* 5^ L^ft); «<*•'); e B (t),0 



with initial conditions 

u (s,s,q,p) =det{e x (s) + e y ) 1/2 
u n (s,s,q,p) = 0, n > 1. 

Corollary 1. Under the additional assumption 

ll( 5 (" g ,p) /l )( < :'7 J P)llL~(RxR^) < oo, 

2<|ct|<4d+7 
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we have the following result on the Ehrenfest timescale T(e) — Ct log(e 1 ): 



sup 

-T{e)<s,t<T(e) 



U £ (t,s)-1 £ (K^;u ;& x (t),e v ^ < C(C* t )£ 1 ~ p(Ct) , 



where p{Ct) can be made arbitrary small, if Ct is chosen small enough. 
Remark 8. 

1. We recall that Z(t,s, q, p) = (i (G^)^ 1 id)^*' 3 ' {q,p)^(-iQ x (t) id) f , thus 
the dependence of Z(t, s, q,p) on q and p is only via the Jacobian of the 
flow. 

2. The expression for the leading order symbol is 

l -i//ii(r, K <T ' s) ( 9 ,p))dr 

U (t,s,q,p) = {det<d y Z(t,s,q,p)) 2 e - 

where the branch of the square root is defined by continuity in t starting 
from uq(s, s,q,p) — det (<d x (s) + O 1 *) 5 , compare Remark [3] The corre- 
sponding FIO is known as the Herman-Kluk propagator in the chemi- 
cal literature. Notice that the dependence of u on q and p is only via 
F k( ' ' (q,p). Likewise, the (<7,p)-dependence in Uk is only via _F K< ' (q,p) 
and its derivatives with respect to q and p. 

3. As an easy corollary we get that the FIO defined in the last theorem 
approximately inherits the properties of U(t,s), i.e. it is almost unitary 
in the sense that 



I s («;<*••>; tt;0 x (t), 9") 2* ^K it ' s) ;u;Q x {t),Q v Y - id 

i £ (V'-^^e^e^V U ( ' t ' s) ;u;e x {t),e y ^ -id 



<C N e N+1 
<C N e N+ \ 



where u = Y2 n =a e ™ u n an d it almost fulfills the group property, i.e. 



u; id, id) 1 E (k^''\ u; id, id) -I e U (t ' s \ u; id, id) 



< C' N e N+1 . 



The result also holds for general <d x and y . The possibility of stating the 
correct dependence of the symbol on the matrices is left to the reader. 



4. In the case of a linear flow 



the approximation becomes exact as the symbols v n in Propositions 
and [3] vanish for n > 2. This can be rephrased as 

1 s (k^;u ; <3 x , 0") = Met E (q,p)) , 

where Met £ denotes the e-dependcnt mctaplcctic representation. 
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5. The Ehrenfest result does not generalise to the higher order corrections. 
This is in accordance with the situation in |BaGrPa99j, where holomor- 
phicity of the symbol is required to get higher order versions of Egorov's 
Theorem on the Ehrenfest timescale. 

6. The proof will produce the following byproduct: the so-called Frozen 
Gaussian Approximation, which is obtained by still choosing «(*' s ) as the 
Hamiltonian flow but keeping uq constant for all q,p and t, is an asymp- 
totic solution of the Schrodinger equation of order 0(e). Thus it approx- 
imates the unitary group for the short times of order e. It will not be 
a valid approximation for longer times because of the uniqueness of the 
symbol. 

Proof. By Theorem [TJ an FIO associated to a C 1 family k^ 1 ' 3 ' of canonical 
transformation of class B and (x, y)-independent symbol u — Yln=o £ ™ u n' u n £ 
C 1 (M, S[Q; 2d]) leaves <S(R d , C) invariant. Thus, we can plug such an operator as 
an ansatz into the time-dependent Schrodinger equation (JTJ) . By Propositions^ 
and [3] we have a representation 

[iej t - of (ho + eh)) I e f 5^E n u„;e x (t),e^ 

/ N+l \ 

\ n=0 ) 

on 5(R rf ,C). We will show that the v n , < n < N + 1 vanish, if and 
u n G S[0,2d], < n < N are chosen properly. Moreover, it will turn out that 

v £ N+2 is of class S[0;3d]. Thus, by TheoremUl l e (k^; J2n=o e " u «; &x (*)> Qv ) 
is an asymptotic solution of order + 2. The statement will then follow by 
Lemma [TJ 

Combining Propositions [2] and [3l one recognises vq as the product of uq and 

+ J t XK ^^P) ■ 2 K<t ' S W) - (i, K (t ' s) (<z,p))) . (34) 

As we do not expect X s (k; 0; O*, <d v ) — to be a good approximation of U (t, s), 
we require (|34|) to vanish. By combining its derivatives with respect to p and q, 
it is easily seen that this is the case if and only if k^' 3 ' is the Hamiltonian flow 
associated to h . 

There are several parts which contribute to v±: 

1 . the zeroth order terms of Propositions [2] and [3] applied to iti , 

2. the first order terms of Propositions [5] and [3] applied to uq, 

3. the zeroth order term of Proposition [5] applied to uq for the subprincipal 
symbol hi. 
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Thus we get the following expression for v± 

d „,(*,«)/ v d 



div. 



(35) 
(36) 



J t S K(t %, P ) + |* K< " W) ■ S K<t ' S> fop) - h (t, K^(q,p) 
^ (d x h + iQ x (t)d e h ) (t, (q, p)) ] f Z- 1 (t, s, q, p)u } 

+ div 2 ^^a" tt '* ) ( 9j p)-te a (t)^ < *"* 5 (g,p)yz-Hi,a,g,p)«o^ (37) 
- Uo itr (z- 1 ^, s, g ,p)5 2 [(9^o + iQ x (t)d e h ) (t, « (M) (g,p) 
+ ^«o - |«otr fz- 1 (t,s ! g,p)X 2 4 (g,p)— G 35 ^) 



- u ^i f*, K (t ' s) (g,p) 
= - ^"otr f-Z -1 ^, s, q,p)^Z(t, s, q,p)j - u hi (t, K (M) (q,p) 

as (|3"5]l . ([35]) and (|3T|) vanish because of the choice of k^'^ as the Hamilton flow. 
As the linearisation of det(A) is det(A)tr(A~ 1 dA) for invertible A, the equation 
v\ — with initial conditions that recover identity, is solved by 

u n {t,s,q,p) = (det(Q y Z{t, s, q,p)))^ exp (-i J h x (r, K (T ' s) (<?,p)) rfrj , 

which is of class S[0; 2d]. 

v n is given by the following expression: 

i^w„_i - ^w„_itr (z^it, s, q,p)-^Z(t : s, q,p)j - m„-i^i (i, K (M) (<7,p)) 



J]div z na 2 , ^-e K (t)z- 1 (i, s , g ,p) efcU „_ 2 ) ) z-^s^p) 
^ i 3 -[ho(*); e B (t),e* 



Un- 



it— J 



k=l 

n-2 

— i 

n— 3 

- i ^ l,-M*); ^ (M) ; e*(<). ©>»-,•-!, 

where we already dropped the terms analogous to (j31)|) - (|3"T)) . The equation 
v n = is easily solved by variation of the constant and has a solution in S[0; 2d\. 
We finally note that the highest order symbol is of class 5[0; 3d]. Thus, we have 
established that the constructed FIO is an asymptotic solution of order N + 2 
on the class of Schwartz functions, so the result follows by the strategy outlined 
at the beginning of the proof. 

We turn to the uniqueness. Assume that there are k^'^ and u <E 5[0;2ef] such 
that 

U s (t, s) - T u; & x (t), 6") < C'{T)e. 
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In this case we have 



; uo;e B (t),e*) -1 E fa t ' a '*-,u;& c (t),Q v } < (C(T) + C'(T)) e 

and thus we get «;(*> s ) = k^' 3 ' and u — u on suppu = M 2d by Proposition |U 
The uniqueness of the higher order corrections follows inductively by the same 
kind of argument. □ 

Proof, (of Corollary[T]) To extend the result to the Ehrcnfcst timescale, we have 
to study the dependence of the remainder's symbol v\ on T = T(e) and to show 
that the growth of 

|a|<4rf+l 

does not exceed O(e) in the e — > limit if Ct is sufficiently small. 
The dependence comes from the elements of F k( ' (q,p) and its derivatives. By 
(fT4"|) . they allow for a bound of the form C'(CT)e~ p ^ CT \ where p(Ct) can be 
made arbitrary small if Ct is chosen small enough. Moreover, v\ has polyno- 
mial growth in these quantities, which follows from the form of the differential 
operators of Proposition [21 the explicit expression of uq and the bound away 
from zero of the determinant of Z(t, s,q,p), compare the proof of Lemma [2] 
Combining these facts, the result follows. □ 



A Gaussian integrals with non-real matrices 

We consider the convex cone C of complex symmetric matrices with positive 
definite real part. Every matrix of C is invertible with its spectrum included in 
the open half plane {z|5Rz > 0}. It follows from matrix theory (see | JoOkReOl] ) 
that each element of C admits an unique square root in C. Furthermore, the 
square root of M is given by the Dunford- Taylor integral (see Ka66j I. §5. 6) 

M x /2 = J_ f zWm-z)- 1 dz 

where the integration path is a closed contour in the half-plane {z|3?z > 0} 
making a turn around each eigenvalue in the positive direction and the value 
of z 1 / 2 is chosen so that it is positive for real positive z. As a consequence, 
the square root M 1 / 2 is an holomorphic function of M. If one considers the 
computation of the Gaussian integral 

j- \ , /0 / e ~^ x ' x dx, 

it is well-known that its value is given by (detM) 1 / 2 = det(M 1//2 ) for positive 
definite real symmetric M. From the above discussion, it directly follows that 
this property extends to any matrix M £ C (see Appendix A in |Fo89j or Section 
3.4. in |H583j for an alternative explanation). 



2G 



References 

[BaGrPa99] D. Bambusi, S. Graffi and T. Paul, Long time semiclassical approx- 
imation of quantum flows: a proof of the Ehrenfcst time, Asympt. Anal., 
21, pp. 149-160 (1999). 

[BiRoOl] J.-M. Bily and D. Robert, The Semi-classical Van-Vlcck Formula. 
Application to the Aharonov-Bohm Effect, Graffi, Sandro (ed.) et al., 
Long time behaviour of classical and quantum systems, Proceedings of the 
Bologna APTEX international conference, Bologna, Italy, September 13-17, 
1999, Singapore: World Scientific, Ser. Concr. Appl. Math., 1, pp. 89-106 
(2001). 

[Bu02] J. Butler, Global h Fourier integral operators with complex-valued phase 
functions, Bull. London Math. Soc., 34 (4), pp.479-489 (2002). 

[Co68] J.D. Cole, Pertubation Methods in Applied Mathematics, Blaisdell Pub- 
lishing Co., Waltham, Mass. (1968) 

[CoRo97] M. Combescure and D. Robert, Semiclassical spreading of quantum 
wave packets and applications near unstable fixed points of the classical 
flow, Asympt. Anal. 14, pp. 377-404 (1997). 

[Fo89] G.B. Folland, Harmonic Analysis in Phase Space, Annals of Mathematics 
Studies 122, Princeton University Press (1989). 

[Fu79] D. Fujiwara, A construction of the fundamental solution for the Schro- 
dinger equation, J. d'Analyses Math., 35, pp. 41-96 (1979). 

[Ha85] G. Hagedorn, Semiclassical quantum mechanics I : The h — > Limit for 
Coherent States. Commun. Math. Phys., 71, p. 77-93 (1980). 

[Ha98] G. Hagedorn, Raising and Lowering Operators for Semiclassical Wave 
Packets, Ann. Phys., 269 (1), pp. 77-104 (1998). 

[HaRoGr04] C. Harabati, J.M. Rost and F. Grossman, Long-time and unitarity 
properties of semiclassical initial value representations, J. Chem. Phys., 
120 (1), pp. 16-30 (2004). 

[Hc75] E.J. Heller, Time-dependent approach to semiclassical dynamics, J. 
Chem. Phys. 62(4), pp. 1544-1555 (1975). 

[He81] E.J. Heller, Frozen Gaussians: a very simple semiclassical approxima- 
tion, J. Chem. Phys. 75(6), pp. 2923-2931 (1981). 

[HcK184] M.F. Herman and E. Kluk, A semiclassical justification for the use of 
non-spreading wavepackets in dynamics calculations, Chem. Phys., 91 (1), 
pp.27-34 (1984). 

[H683] L. Hormander, The Analysis of Linear Partial Differential Operators I, 
Springer- Verlag, New York (1983). 

[JoOkReOl] C.R. Johnson, Y. Okubo, R. Reams, Uniqueness of Matrix Square 
Roots and an Application, Lin. Alg. Appl, 323, pp. 51-60 (2001). 



27 



[Ka94] K. Kay, Integral expressions for the semi-classical time-dependent prop- 
agator, J. Chem. Phys., 100 (6), pp.4377-4392 (1994). 

[Ka06] K. Kay, The Herman-Kluk approximation: Derivation and semiclassical 
corrections, Chem. Phys., 322, pp.3-12 (2006). 

[Ka66] T. Kato, Perturbation Theory for Linear Operators, Springer- Verlag, 
New York (1966). 

[Ki82] H. Kitada, A calculus of Fourier integral operators and the global funda- 
mental solution for a Schrodinger equation, Osaka J. Math., 19, pp. 863-900 
(1982). 

[KiKu81] H. Kitada and H. Kumano-Go, A family of Fourier Integral Operators 
and the fundamental solution for a Schrodinger equation, Osaka J. Math., 
18, pp.291-360 (1981). 

[LaSiOO] A. Laptev and I.M. Sigal, Global Fourier Integral Operators and semi- 
classical asymptotics, Review of Math. Phys., 12 (5), pp. 749-766 (2000). 

[Ma02] A. Martinez, An Introduction to Semiclassical and Microlocal Analysis, 
Universitext, Springer- Verlag, New York (2002). 

[RoSw07] V. Rousse, T. Swart, Global L 2 -Boundedness Theorems for Semiclas- 
sical Fourier Integral Operators with Complex Phase, submitted, preprint 
is available online at http://arxiv.org/abs/0710.4200. 

[Ta04] D. Tataru, Phase space transforms and microlocal analysis, in Phase 
space analysis of partial differential equations Vol. II (Pubbl. Cent. Ric. 
Mat. Ennio Giorgi, Scuola Norm. Sup. Pisa), pp. 505-524 (2004). 



28 



